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ABSTRACT 



This study was conducted to establish the performance parameters 
of a radial inflow, dual discharge turbine and to determine the effect 
of variations in axial clearance on these parameters . The representa- 
tive stream surface is taken at the outer discharge radius instead of 
at a computed mass-average discharge radius, as was done previously. 

This technique results in considerably simplified computations and in 
better correlation of the rotor loss parameters. 

Tests were conducted at axial clearances from 0.015 to 0.081 inches 
and at total-to-static pressure ratios from 1.2 to 1.7 for each clear- 
ance. The test installation is located at the Turbo-Propulsion Labora- 
tory of the Naval Postgraduate School, Monterey, California. 
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TABLE OF SYMBOLS 



Symbol 


FORTRAN 


Definition 


Units 


A i 


- - - - 


Meridional cross-sectional 
area of rotor inlet 


ft 2 


A 5 


- - - - 


Cross-sectional area of five 
inch inlet pipe 


ft 2 


B 1 


B1 


Perpendicular distance between 
shroud inner extremities 


in 


C 


— “ “ “ 


Flow measurement factor dependent 
on units, orifice diameter, and 
type of pressure taps 


- - - 


C fl 


CF1 


Pressure conversion factor 


lb/ft 2 /in Hg 


Cine 


CLNC 


Rotor axial tip clearance, 
on one side of blade only 


in 


c o 


CO 


Velocity corresponding to the 
isentropic enthalpy drop 
through the turbine 


ft/sec 


c 

p 


CP5 


Constant pressure specific heat, 
based on total inlet temperature 


BTU/lb -°F 
m 


Cp(avg) 


CP 


Constant pressure specific heat, 
based on average temperature t 


BTU/lb -°F 
m 


D 1 


D1 


Rotor diameter at inlet, 9.40 in. 


in 


°20 


- - - - 


Rotor outer discharge diameter, 
5.88 in. 


in 


F 


SR 


Scale reading 


lb 


Fr 


- - - - 


Reynolds number correction 
factor 


- - - - 


°Hg 


GHG 


Specific gravity of mercury 
at t 

rm 


- - - - 


G H20 


GWR 


Specific gravity of water 
at t 

rm 


— — — — 


g 


32.174 


Unit Conversion Factor 


ft-lb /lb-sec 2 
m 


h e 


----- 


Total enthalpy at equivalent 
state point E 


BTU/lb 

m 
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Symbol FORTRAN 



Hp f 


HPF 


h r 


- - - 


H refl 


HREF1 


h l 


- - - 


CM 

-C 


- - - 


(avg) 


HI 


AH. 

is 


DHIS 


AH 

w 


- - - 


Ah, 

lvc 


DPVC 


J 


778.16 


k. 

is 


HEAD 


M 


TNET 


M f 


BFM 


M ref 


TNETR 


M 1 


ACH1 


N 


RPM 


N * 
ref 


RPMR 


P 

atm 


PATM 


p t0 


PTO 



Definition 

Power required to overcome 
bearing friction 

Relative total enthalpy at 
rotor inlet 

Measured reference pressure 
applied to water manometer 
board 

Static enthalpy at rotor 
inlet 

Static enthalpy at rotor dis- 
charge 

Average static pressure ahead 
of rotor 

Isentropic total enthalpy drop 
across turbine 

Enthalpy drop representing 
work output of turbine 

Measured pressure differential 
across flow orifice, from vena 
contracta taps 

Joule's constant 

Isentropic head coefficient 

Net torque output of turbine 

Bearing friction moment 

Referred net torque output 
of turbine 

Mach number at rotor inlet 

Measured turbine speed 

Referred turbine speed 

Barometric pressure 

Absolute total pressure ahead 
of turbine 



Units 

h P 



BTU/lb 

m 



in Hg gage 



BTU/lb 

m 

BTU/lb 

m 

in H 2 0 gage 
BTU/lb 

m 

BTU/lb 

m 



in Hg 



ft-lb/BTU 

ft-lb 

ft-lb 

ft-lb 



rpm 
rpm 
in Hg 
in Hg abs 
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Symbol 


FORTRAN 


Definition 


Units 


p t0 /p l 


PR1 


Pressure ratio across scroll 
and guide vanes 


- - - - 


P t0 /p 2 


PR2 


Pressure ratio across turbine 


- - - - 


P E 


- - - - 


Equivalent rotor inlet static 
pressure 


in Hg abs 


P 0 


PS 5 


Absolute static pressure at 
turbine inlet 


in Hg abs 


P 1 


PI 


Average absolute static 
pressure at rotor inlet 


in Hg abs 


p ivc' 


PUVC 


Measured pressure upstream of 
flow orifice , from vena 
contracta tap 


in Hg gage 


P 1VC 


PVC 


Absolute pressure upstream of flow 
orifice from vena contracta taps 


in Hg abs 


P 2 


P2 


Static discharge pressure 


in Hg abs 


p 5 ' 


P5P 


Measured static pressure at 
turbine inlet 


in Hg gage 


R 

g 


53.3448 


Gas constant for air 


ft-lb/lb - °R 
m 


R i 


Rl 


Rotor inlet radius, 4.70 in. 


in 


R 20 


R20 


Rotor outer discharge radius, 
2.94 in 


in 


r 


DR 


Isentropic degree of reaction 


- - - - 


T 


T 


Torque measured by dynamometer 


ft-lb 


TCD 


TCD 


Torque calibration data 


counts 


TQ 


TQ 


Torque indicator reading 


counts 


t e 


- - - - 


Equivalent rotor inlet 
temperature 


°R 


T ref 


TR 


Referred torque output for 
no-loss case 


ft-lb 


T to 


T5 


Total temperature ahead of 
turbine 


°R 


T o 


TO 


Static temperature ahead of 
turbine 


°R 
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Symbol 


FORTRAN 


Definition 


Units 


T i 


Tl 


Static temperature ahead of 
rotor 


°R 


V 




Static temperature at rotor 

inlet for isentropic expansion 

from P - to p, 
to 1 


°R 


T * 
1 




Static temperature at rotor 
inlet resulting from rotor 
incidence losses 


°R 


T ' 
2 


“ “ “ “ 


Static discharge temperature 

resulting from isentropic 

expansion from p to p 
h 2 


°R 


rp f I 

2 


***** 


Static temperature at discharge 
resulting from isentropic 
expansion from P ^ to p^ 


°R 


T 4 


T4 


Total temperature ahead of 
flow orifice 


°R 


T 20 


T20 


Lubricating oil inlet 
temperature 


°F 


T 

21 


T21 


Lubricating oil temperature 
at discharge from right bearing 


°F 


T 

22 


T22 


Lubricating oil temperature 
at discharge from left bearing 


°F 


A T. 
is 


DTIS 


Isentropic temperature drop 

from P A to p. 
tO v 2 


°R 


AT w 


DTW 


Temperature drop corresponding 
to actual work output of turbine 


°R 


t 


A 


Average temperature through 
turbine 


°F 


tare 


STARE 


Tare of precision scales 


lb 


t 

cj 


TCJ 


Thermocouple reference cold 
junction temperature 


°F 


t 

rm 


TRM 


Control room temperature 


°F 


u i 


Ul 


Peripheral speed of rotor 
at inlet 


ft/sec 


U 20 


U20 


Peripheral speed of rotor 
at outer discharge radius 


ft/sec 
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Symbol 


FORTRAN 


Definition 


Units 


V 


- - - " 


Thermocouple reading 


mv 


v o 


VO 


Velocity at turbine inlet 


ft/sec 


V 1 


VI 


Absolute flow velocity at 
rotor inlet 


ft/sec 


V ul 


VUL 


Peripheral component of 


ft/sec 


V ml 


VM1 


Meridional component of 


ft/sec 


V ith 


VITH 


Theoretical absolute rotor 
inlet velocity for no losses 


ft/sec 


V 20 


V20 


Absolute discharge flow 
velocity at outer radius 
of rotor 


ft/sec 


V u20 


VU20 


Peripheral component of V^q 


ft/sec 


V m20 


VM20 


Meridional component of 


ft/sec 


V 4 


V4 


Thermocouple reading ahead 
of flow orifice 


mv 


V 5 


V5 


Thermocouple reading at 
turbine inlet 


mv 


o 

CM 

> 


V20 


Thermocouple reading at 
lube oil inlet 


mv 


V 21 


V21 


Thermocouple reading at 
lube oil discharge from 
right bearing 


mv 


V 22 


V22 


Thermocouple reading at 
lube oil discharge from 
left bearing 


mv 


W 1 


W1 


Relative flow velocity at 
rotor inlet 


ft/sec 


W ul 


WU1 


Peripheral component of 


ft/sec 


W 20 


W20 


Relative flow velocity at 
rotor discharge outer radius 


ft/sec 


W u20 


WU20 


Peripheral component of 


ft/sec 


W 20th 


W20TH 


Theoretical relative flow 
velocity at outer discharge 
radius, for isentropic expansion 
from p^ to p^ 


ft/sec 
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Symbol 


FORTRAN 


Definition 


Units 


4 


WVC 


Mass flow rate* for vena 
contracta taps 


lb /sec 
m 


“ref 


WVCR 


Referred mass flow rate 


lb /sec 
m 


w 

oil 


- - - - 


Mass flow rate of lubricating 
oil 


lb /sec 
m 


Y 1 


Y 


Factor accounting for 
compressibility effects in 
flow measuring orifice 


- - - - 
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Units 



Symbol FORTRAN Definition 

GREEK SYMBOLS 

CK A Area multiplier accounting for 

thermal expansion of flow orifice 



C* 1 


ALPH1 


Absolute rotor inlet flow 
angle 


degrees 


o 

CM 


ALF2 


Absolute rotor discharge 
flow angle at outer radius, 
for no bearing loss case 


degrees 


° < 20L 


ALF2L 


Absolute rotor discharge flow 
angle at outer radius, bearing 
losses considered 


degrees 


Pi 


BETA1 


Relative rotor inlet flow 
angle 


degrees 


o 

CM 

(XL 


BET20 


Relative rotor discharge 
flow angle at outer radius 


degrees 


Y 


GAMMA 


Ratio of specific heats 

based on T . 

to 


- - - - 


^(avg) 


GAM 


Ratio of specific heats 
based on t 


- - - - 


(Y-D/Y 


EXP 


Exponent 


- - - - 


% 


DEL 


Referred pressure ratio, 

/1A.7 

to 


- - - - 


£ 


EPS 


Factor which corrects flow 
rate through the turbine 
for varying values of V 




n 


ETA 


To tal-to- static turbine 
efficiency, bearing losses 
not considered 


per cent 


% 


ETAL 


Total- to-sta tic turbine 
efficiency, including 
bearing losses 


per cent 


d 


PHI 


velocity coefficient for 
scroll and guide vanes 


- - - - 


V 


PSI 


Velocity coefficient for 
rotor, based on ratio of 
actual to theoretical 
relative discharge velocities, 
bearing losses not considered 
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Symbol 


FORTRAN 


Definition 


Units 


Yl 


PSIL 


Velocity coefficient for 
rotor, bearing losses 
included 


- - - - 


fo 


RHO 


Density at turbine inlet 


lb /ft 3 
m 


fl 


RHO 


Density at rotor inlet 


ib /ft 3 
m 


© 


THETA 


Referred temperature ratio, 

T /518.7 
to 


- - - - 


lj 


ROSP 


Rotational speed of rotor 


radians/sec 


(PN 


ZETAN 


Loss coefficient for scroll 
and guide vanes 


- - - - 


ipR 


ZETR 


Rotor loss coefficient, 
bearing losses not con- 
sidered 


- - - - 


|pRL 


ZETRL 


Rotor loss coefficient, 
with bearing losses 


- - - - 
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SECTION I 



INTRODUCTION 

As requirements for smaller more powerful turbines increase, to keep 
pace with the needs of the rapidly expanding aerospace industry, more and 
more attention is given to radial turbines • Since the radial turbine can 
extract more energy per stage than an axial turbine, it is more useful in 
applications where space limitations preclude the use of a staged axial 
turbine. Currently, it is therefore used extensively for aircraft elec- 
tric power generation, in aircraft cockpit and equipment cooling systems, 
in cryogenic systems, and in missile auxiliary power systems. 

As a consequence of the increasingly widespread use of the radial 
turbine, the establishment of performance parameters as an aid to the 
design of these turbines is of fundamental importance. It is the objec- 
tive of this report to establish some of the performance parameters for 
a particular radial inflow, dual discharge turbine using air as an 
operating fluid and relating these parameters in such a way that they may 
be applied to geometrically similar turbines operating with different 
fluids and in different environments. 

Although similar studies have previously been conducted on this 
machine, this investigation differs in two important aspects. First, 
several significant and necessary changes in the internal dimensions of 
the turbine were made, thus resulting in a somewhat different machine 
than that previously tested. Also, the method of analysis used is a 
considerably simplified one, which yields more significant results with 
less computational effort than those methods previously employed by re- 
ferring all losses to the outer discharge radius rather than a computed 
average radius. 
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Additionally* this project had as an objective a significant expan- 
sion of the operating ranges examined beyond those previously considered. 

Because the evaluation of turbine performance parameters involves a 
formidable computation requirement* the data obtained in the experiments 
were reduced by means of the computer programs described herein, which are 
written in FORTRAN IV language for the IBM 360 computer system used at 
the Naval Postgraduate School. 

The author gratefully acknowledges the guidance and counsel of 
Professor Robert D. Zucker* of the Department of Aeronautics, Naval 
Postgraduate School, and of Professor Michael H. Vavra of the same 
department. Thanks are due also to Mr. James E. Hammer and his coworkers 
on the technical staff of the Turbopropulsion Laboratory for their many 
hours of hard work in modifying, repairing, and operating the turbine 
test rig. 



18 



SECTION 2 



THE loC.Po RADIAL TURBINE 

The turbine used in this investigation is a dual discharge, radial 
inflow turbine installed at the Turbopropulsion Laboratory of the Depart- 
ment of Aeronautics, Naval Postgraduate School* It was originally de- 
signed and built by the Aerojet General Corporation, for application as 
the power supply for a missile propulsion combustion cycle* The origi- 
nal performance evaluation was conducted by Vavra (1), with subsequent 
investigations by Finn (2), Riley (3), and Boshoven (4)* 

2*1 Description of Turbine Assembly 

The basic turbine test installations used are shown in Figures 1 
and 2* The turbine assembly consists of % the wooden inlet casing, with 
its inner contour in the shape of a varying diameter torus, or scroll; a 
circular stainless steel guide vane assembly consisting of seven circular- 
arc blades and the rings themselves; the rotor assembly; and two plexi- 
glass shrouds* 

The inner contour of the inlet casing was formed by casing plaster 
of paris around a wooden insert of the desired scroll dimensions* The 
inner surface of the scroll was varnished to prevent erosion. The wooden 
insert used in this process is shown in Figure 3, and Figure 4 presents 
a cross-section of the wooden casing, showing the shape of the scroll and 
the locations of the various scroll pressure taps* It is worthy of note 
that the flow entry is not tangential, as is customarily the case in 
radial turbine installations* This characteristic is a result of design 
requirements in the originally proposed missile installation* 
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A cross-section of the turbine is shown in Figure 5. In this view 
are indicated the turbine dimensions and minimum clearance as they 
existed for the test runs. The clearance , or axial distance between the 
rotor and the shrouds* was increased during the tests by means of circu- 
lar metal shims inserted between the shroud flanges and the aluminum 
rings on which they rest. 

The series of tests conducted utilized two different rotors. The 
"dummy" rotor, so called because it does not rotate, is a special device 
designed to permit measurement of the torque generated by the incoming 
flow and thus provides a measure of scroll and guide vane performance. 

This rotor is shown in Figure 6. It has a diameter of 9.50 inches, an 
axial length of 8.50 inches, and 36 meridional blades. This diameter 
places the blade leading edges at the radius which corresponds to that of 
the annulus containing the static pressure taps ahead of the rotor, and 
the axial length extends the blade trailing edges beyond the shrouds. 

The "active" rotor, or the power-producing rotating element, is 
shown in Figure 7. It is 9.40 inches in diameter and has an axial length 
of 5.00 inches. The rotor discharge radius is 1.76 inches at the hub and 
2.94 inches at the tip. The fifteen blades have a spacing at the dis- 
charge varying from 0.737 to 1.232 inches, hub to tip respectively. The 
actual rotor discharge area on one side of the rotor is 6.43 square inches, 
including an effective flow area through the radial gap between the blades 
and the shroud taken at the discharge plane (1). A view of the left 
rotor discharge, showing the plexiglass shroud and several of the eight 
shroud pressure taps which measure static pressure at the rotor inlet, 
is given in Figure 8. 
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2.2 Modifications to Turbine Assembly 

When the turbine assembly was dismantled in preparation for this 
series of tests * it was discovered that one of the plexiglass shrouds 
was deeply gouged. Since removal of enough material to eliminate this 
fault would have destroyed the contour of the shroud * it was decided to 
cut the contours of both shrouds so”that they would exactly fit the rotor 
contour. Additionally* it was decided to correct several dimensional 
discrepancies noted in earlier reports. 

Riley (3) indicates a distance of 0.943 inch between the inner ex- 
tremities of the shrouds. Since the axial width of the blade tips is 
0.894 inch* this indicates a basic (unshimmed) axial tip clearance of 
0.025 inch on each side of the blade. Roshoven (4)* however* determined 
the distance between shrouds to be 0.966 inch* with a resulting unshimmed 
clearance of 0.036 inch. Additionally* Eoshoven discovered that the in- 
side extremity of the left shroud (as viewed from the turbine inlet) was 
offset 0.016 inch more from the guide vane ring than the right shroud. 

In order to correct these discrepancies* an extensive series of 
measurements was undertaken. The guide vane rings* which were found to 
be warped* were lapped until they were flat* and all rough edges* burrs, 
and high spots on the inner surfaces of the casing halves were carefully 
removed. This resulted in a consistency to within 0.0025 inch in the 
dimension between the shroud support rings* at all points around the pe- 
riphery. Initial measurements showed a maximum variation of 0.01 inch 
in this dimension. 

Measurements with the reworked shrouds installed then showed a 
minimum clearance of 0.0477 inch* with the left shroud offset axially 
outward 0.0056 inch more than the right. The shroud flange thicknesses 
were then reduced so as to compensate for this dissymmetry and to yield 
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a final unshimmed clearance of 0.015 inch. Subsequent measurement of the 
axial travel of the active rotor*, from contact at one shroud face to con- 
tact at the other*, verified that all calculated clearances used in these 
tests were correct to within + 0.002 inch. 

Additional modifications were made to the dummy rotor in order to 
correct the axial stagger and unequal widths of the blades and to correct 
as much as possible the large gaps which were found to exist between the 
rotor blade contours and the shroud contours. By careful machining, the 
blade contours were made uniform*, and a much closer fit was made to the 
shroud contours*, even though it was not possible to completely correct 
this condition and still maintain the blade width. 

It is thus apparent that the results obtained from the turbine after 
this series of modifications could differ to a marked degree from those 
previously obtained under identical operating conditions. It also seems 
to this author that both Riley (3) and Boshoven (4) were operating at 
clearances somewhat greater than stated in their respective studies. 

Based only on these modifications, the distance between the shroud inner 
extremities for their tests would have been approximately 0.989 inch, 
yielding a minimum clearance of 0.048 inch. This may be compared with 
the minimum clearance of 0.0245 inch reported by Riley and 0.036 inch 
reported by Boshoven. 

Other modifications to the components of the test rig were made and 
are covered in the sections to which they apply. In this section only 
those modifications basic to the radial turbine itself are discussed. 
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SECTION 3 



ANALYSIS OF SCROLL AND GUIDE VANE PERFORMANCE 

In order to evaluate the performance parameters of a turbine, it 
is necessary to divide the losses into those which occur in the inlet 

guide vanes and those which occur in the rotor itself* Riley (3) and 
Boshoven (4) have used the special dummy rotor installation previously 
described (Figure 2) to determine the losses in the stator, or scroll 
and guide vane assembly* The installation determines the average condi- 
tions at the rotor inlet and thus provides an indication of scroll and 
guide vane efficiency as well as the absolute velocity (in both mangi- 
tude and direction) at the rotor inlet* 

Riley (3) concluded that the velocity coefficient, which is a measure 
of the losses occurring in the inlet manifold and guide vanes, as well as 
the absolute inlet flow angle, was essentially constant in the range of 
turbine pressure ratios and axial clearances examined. The results ob- 
tained by Boshoven (4) indicate that this assumption is invalid; in 
fact, the velocity coefficients were found to increase with pressure ratio 
for a given clearance, and the absolute rotor inlet flow angles were found 
to decrease linearly with pressure ratio* Specific results of these ear- 
lier tests are discussed in detail in Section 3*3* 

It was originally intended to utilize the dummy rotor test data of 
Boshoven (4) in this experiment. However, because of the extensive modi- 
fications previously discussed, it was necessary to repeat these tests. 
This repetition yielded the additional advantage of permitting testing 
over a wider range of clearances* 
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3.1 Instrumentation and Test Procedures 



The dummy rotor shown In Figure 6 is mounted on self-aligning ball 
bearings and replaces the turbine rotor. Its purpose is to turn the flow 
from the guide vanes Into the axial direction and to discharge it with no 
peripheral velocity component. If this condition is met, the losses in 
the scroll and guide vanes may be determined by measuring the torque ex- 
erted on the static dummy rotor, the flow rate, the conditions at the 
manifold inlet, and the static pressure at the dummy rotor inlet. 

As was previously mentioned, the dummy rotor has an axial length of 
8.50 inches, which extends the blade trailing edges beyond the shrouds. 

To Insure that the discharge flow had no whirl components, flow etraight- 
eners, also shown in Figure 6, were used. These devices consist of alumi- 
num caps containing a 1.00 inch depth of 0.125 by 0.0015 inch honeycomb 
material. The caps slide over the rotor shaft and cover the rotor blade 
trailing edges. The entire dummy rotor assembly is supported by brackets 
attached to the manifold casing, as shown by Figure 2. 

Air for turbine operation was supplied by an Allls-Chalmers 12-stage 
axial compressor. Air from the main supply line passes through a four- 
inch pipe, a settling tank, a five-inch pipe, and then into the turbine. 
The flow rate was regulated by remotely controlled butterfly valves in the 
inlet line to the turbine and m the outside discharge line of the Allls- 
Chalmers compressor. (See Figure The controls for operating these 
valves are located in the control room (Figure 10). 

The flow rate of the incoming air was measured with a sharp-edged 
orifice located in the four-inch inlet pipe. This orifice is also shewn 
in Figure 9. The pressure ahead of the orifice and the pressure differ- 
ential across the orifice were measured with standard vena eontraeta 
taps. The total temperature ahead of the orifice was measured with an 



iron-cons tantan thermocouple installed in a Kiel temperature probe. The 
orifice installation conforms to standards set forth by Stearns, et^ al . 
(5). 

The total inlet conditions in the five-inch pipe just upstream of 
the turbine were measured with a Kiel probe * also incorporating an iron- 
constantan thermocouple . The static pressure at the same location was 
measured by two static pressure taps in the line. The static pressure 
at the rotor inlet was measured by pressure taps located at eight pe- 
ripheral stations in each shroud. The axial location of these taps is 
depicted in Figure 5* and their peripheral locations are shown in Figure 
4. These taps were moved to the positions shown from their original loca- 
tions near the shroud inner extremities by Boshoven (4) in order to re- 
duce the effects of local pressure perturbations experienced in previous 
tests. 

The torque exerted on the dummy rotor by the flow was measured on a 
25-pound capacity Toledo precision scale „ shown in Figure 2. A twelve- 
inch lever arm was attached to the rotor shaft and incorporated a simple 
provision for leveling* so that the force was exerted on the scale through 
an exactly vertical link. This was done to preserve the direct relation- 
ship between scale reading and moment. A counter balance was added to the 
lever arm in order to extend the capacity of the scales to accommodate 
higher mass flow rates. 

A general schematic of the instrumentation used in this and the sub- 
sequent active rotor tests is shown in Figure 11. From this schematic it 
may be seen that the pressure upstream of the flow measuring orifice, 
the total and static pressures at the turbine inlet* and the reference 
pressures 1 and 2 were read from an eleven column* 40-inch mercury mano- 
meter board located in the control room. This board was referenced to 
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atmospheric pressure* Additionally , one arbitrarily selected scroll pres- 
sure tap was connected to the mercury board, and this pressure, together 
with the total inlet pressure, was used to set the pressure ratio across 
the scroll and guide vanes. The total inlet pressure displayed on the 
mercury board was used only for setting this pressure ratio since a more 
representative total pressure was computed in the data reduction. 

The sixteen shroud pressure taps that measure the static pressure at 
the rotor inlet were read on a fifty-column, 96-inch water manometer board. 
Also displayed on this board were the pressures at various stations in the 
scroll itself, although these pressures were not used in the tests. These 

i 

pressures were measured against atmospheric pressure for the entire series 
of dummy rotor runs. However, additional, controlled, pressures could have 
been imposed on the water board reservoirs by means of regulated house air 
if this had become necessary in order to avoid exceeding the limits of the 
manometer board. These pressures (the reference pressures 1 and 2 men- 
tioned above) were displayed on the mercury manometer board. The first of 
these pressures, displayed on column 10 of the mercury board, was recorded 
as a test datum since it applied to the sixteen columns which indicated 
the rotor inlet pressures. The value of this pressure was thus zero for 
all dummy rotor runs. 

The pressure differential across the flow orifice was measured on a 20- 
inch Meriam model 34FB2 mercury micro-manometer, which has a measuring ac- 
curacy of j^).001 inch. A pair of these instruments is shown in Figure 12. 

The potentials generated by the iron-constantan thermocouples were 
measured on a 48-channel Brown potentiometer, manufactured by Minneapolis- 
Honeywell. This potentiometer, shown in Figure 13, was located in the 
control room and used an ice bath in the test cell as a cold junction 
reference . 
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A series of eighteen test runs were made on the dummy rotor instal- 
lation, of which eight are documented in this report* These eight runs 
represent clearances of 0*015, 0*024, 0.030, 0 * 041 , 0*051, 0 * 061 , 0 * 072 , 
and 0*081 inches, respectively* For each clearance ten pressure ratios 
were sampled, ranging in value from 1.05 to 1*65. These pressure ratios 
represent the ratio of total turbine inlet pressure to the static pres- 
sure at the rotor inlet* 

This large number of runs was necessitated by several factors* Data 
from the initial runs possessed an unacceptable degree of scatter. This 
was attributed to an unstable supply of source air from the Allis-Chalmers 
compressor* When this unsatisfactory trend was discovered, the speed 
control system was overhauled, and the subsequent test results were signi- 
ficantly improved. 

Additionally, the stability of the system was improved by judicious 
use of the aforementioned butterfly control valve. Experimentation showed 
that if less of the excess air from the compressor were discharged to 
the outside, the compressor would be operating at a higher pressure ratio,/ 
(with a greater pressure drop across the valve feeding the turbine), and 
the pressurized air in the house plumbing would act as an additional plenum. 
The system fluctuations would then be substantially reduced* These im- 
provements in equipment and operating technique, as well as the inevita- 
ble improvement in data taking technique, led to such improved results that 
many runs were repeated in order to obtain the best possible test data. 

As first observed by Boshoven (4), it was noted that the pressure 
distribution around the shroud was very sensitive to the relative posi- 
tions of the rotor blade leading edges. To eliminate this factor as a 
possible variable, it was arranged to permit the return of the rotor to 
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the same location each time a new pressure ratio was set. This was ac- 
complished by means of a wing nut on the end of the rod connecting the 
moment arm to the scale. 

It was also observed that the pressures seemed to be higher by about 
4 to 5 inches of water at the right side of the rotor inlet than at the 
left. This condition could be eliminated by moving the rotor approximately 
0.001 inch axially to the left; but since the condition was also present in 
the scroll itself, it was attributed to scroll pressure distribution rather 
than to misalignment. The rotor was centered for all runs. 

3.2 Theoretical Analysis 

From the law of moment of momentum (6), the moment M exerted on the 
dummy rotor by the flow is given by 

M - ( W/g) Ri Vui <« 

where W is the mass flow rate, is the inlet radius of the dummy rotor, 
and is the peripheral velocity component at R^. Equation (1) is 
valid for the case of steady flow in the dummy rotor because the pressures 
acting over the inlet and discharge areas cannot produce a moment about 
the axis, and because the flow straighteners insure an axial discharge 
velocity. 

A measure of the performance of the scroll and guide vanes is given 
in the form of the velocity coefficient, (j) , which is defined by 

<j> = V Vui, < 2 > 

where is the average velocity at the dummy rotor inlet and 

the theoretical velocity at the dummy rotor inlet which would result from 

an isentropic expansion from the manifold inlet total conditions P^q, T^q, 
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to the static pressure p. at the dummy rotor inlet* This process is 

1 

shown on the T-s diagram of Figure 14. The component velocities may be 
seen on Figure 15. 

From the energy equation for an adiabatic , no work process, the 
velocity at the dummy rotor inlet is 

Vi - / 2-g T c P (Tto - Ti ) (3) 

If use is made of the isentropic pressure»temperature relation for a 
perfect gas and of Equation (2), may be alternately expressed as 

Vi = i v /2.gTc p T tl) [l-(-p 1 /fto) J T] ' < 4 > 

Application of the equation of continuity at the dummy rotor inlet 
yields the following expression for the average meridional velocity 
component i 

Vnni = W/ fi Al (5) 

where ^ is the density at the dummy rotor inlet, and is the cross- 
sectional area normal to the flow (meridional cross-section). The area 
A^ is given by 

Ai - Ztt Ri Bi < 6 > 

where is the distance between the shrouds at the rotor inlet. Using 
the equation of state for a perfect gas, ^ may be expressed as 

= T’l/RgTi (7) 

Substitution of this relation and Equation (6) into Equation (5) yields 
the resultant expression for : 
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From Equation (1), the peripheral velocity may be expressed as 




Rj. W 



(9) 



and may additionally be given by 

Vi = m + < 10 ^ 



Thus, three equations expressing the velocity are established, 
these being Equations (3), (4), and (10). Upon substitution of Equations 
(8) and (9) into Equation (10), it may be seen that this set of three 
equations contains three unknowns. The unknown quantities to be deter- 
mined are the velocity V^, the velocity coefficient (j> , and the tempera- 
ture T^, all other quantities being either constants or determined direct- 
ly from measured data. Solution of this system of equations establishes 
the velocity components, which may then be used to calculate the absolute 
flow angle at the dummy rotor inlet, according to the relation 



^ i - tan 1 (Vuj.y/ Vrvn ) 

A nozzle loss coefficient, may be obtained from the relation 

IfN = L - ^ - 1- Tf\ (12) 

where is an efficiency based on the ratio of the actual-to-isentropic 
temperature drop across the scroll and guide vanes. It is introduced here 
only to illustrate the significance of L and 0. 
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3.3 Reduction of Data; Program SCROLL 

The computer program SCROLL, listed in Appendix A, determines the 
unknown quantities, j> , and , and the absolute inlet angle, 

from the test data using an iterative procedure. The program makes 
temperature corrections to the physical properties of air and to the 
specific gravities of the water and mercury used as manometer fluids. 

It converts the manometer readings and thermocouple voltages into pres- 
sures and temperatures, as required, and then establishes the mass flow 
rate. The average conditions at the dummy rotor inlet are then calcu- 
lated from the measured torque in accordance with the previously given 
relations. 

Program SCROLL has been described in detail by Riley (3) and Boshoven 
(4). Except for conversion to the FORTRAN IV language, the program used 
in this experiment remains essentially the same; therefore, the details 
of the computational processes will not be repeated. However, certain 
minor changes were incorporated in the temperature corrections applied to 
the manometer fluids, in the thermocouple equations used, and in the pres- 
sure and flow rate solutions . These changes are briefly described in 
Appendix A, although they become evident if the program used in this re- 
port is compared with either of the two above-listed references. 

3.4 Discussion of Results 

The results of program SCROLL consist primarily of values of the 
velocity coefficient, (j> , and the absolute flow angle at the rotor in- 
let, c These values are plotted versus M^, the Mach number at the 

rotor inlet, as opposed to the pressure ratio, P /p , because the Mach 

to I 

number is a more significant parameter in the application of test results 
to geometrically similar machines. Figure 16 shows as a function of 
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for each of the eight clearances tested, and Figure 17 is a similar 
plot for as a function of M^. Each of the curves was fitted to the 
data obtained for that run by the least squares method. The tabulated 
results for each of the eight runs are given in Appendix A. 

It may be seen from Figure 16 that the velocity coefficient 
increases with Mach number, the average increase being approximately 4.75 
per cent over a Mach number range of 0.2 to about 0.8. This percentage 
increase is practically the same regardless of clearance. This increase 
is due, at least in part, to the decreasing coefficient of friction in 
the scroll and guide vanes, which results from the increasing Reynolds 
number. 

The coefficient <j generally increases with decreasing clearance for 
a given pressure ratio with the maximum values of cf> occurring at the mini- 
mum clearance of 0.015 inches. As clearance is increased beyond a value 
of about 0.070 inches, the incremental drop in <f) per increment of clear- 
ance increase appears to become much larger. An explanation for this oc- 
currence is that as clearance is increased, an increasing amount of flow 
escapes around the blade tips. This flow loses angular momentum by wall 
friction and subsequently strikes the dummy rotor blades at some lesser 
radius, thus imparting a reduced amount of torque to the rotor. The cumu- 
lative effect of this energy loss would logically become more pronounced 
at increased clearances. The overlapping of certain of the curves in Fig- 
ure 16 at the lower Mach numbers occurred because of the relative in- 
sensitivity of the instrumentation at low pressures. 

Figure 17 shows that the absolute rotor inlet flow angle, o<^, in- 
creases with increasing Mach number up to a Mach number of approximately 
0.4. Beyond that point, o<^ decreases with increasing Mach number. These 
changes are small, although they appear to be significant because of the 
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highly expanded scale of Figure 17. In fact* the average maximum change 
in c<^ over the range of Mach numbers is approximately 0.52 degrees. In 
general* increasing the mass flow rate will increase and* of course* 
the components and V The variations in the magnitudes of these 
two quantities as the Mach number is increased determine the variations 
of as shown by Equation (11). 

It is also shown on Figure 17 that the angle increases as clear- 

ance increases. This is primarily because V „ is inversely proportional 

ml 

to the meridional cross section* which increases with clearance since it 
depends on the distance between the shrouds. By Equation (11), the smaller 
values of V ^ result in larger inlet flow angles. 

It was concluded by Riley (3) that the velocity coefficient <f> was es- 
sentially independent of the pressure ratio and that a constant 

representative value of 0.889 was sufficiently accurate. Although he 
determined that o ( ^ decreases with increasing pressure ration* he also 
assumed a representative value of 80.0 degrees for 

Boshoven (4) subsequently investigated the scroll and guide vane per- 
formance in greater detail and obtained results that agree in a qualita- 
tive sense with those described above. Because of the many modifications 
made to the machine components, both in preparation for this experiment 
and during the course of Boshoven 5 s* a more precise comparison is not 
warranted. It is worthy of note* however* that Boshoven did not obtain 
the continuous decrease in ^ with increasing clearance which is apparent 
from Figure 16. Further* the increase in ^ with Mach number found in the 
present tests is more pronounced than Boshoven 9 s, his increase being on 
the order of 2 per cent for a given clearance* over the range of Mach 
numbers. It is felt that these changes have occurred because of the 
machining of the rotor so that it more closely fits the shroud contours* 
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thus minimizing flow loss around the blades and permitting closer mini- 
mum clearances. 

An examination of the data points given in Figure 16 for 0 as a 
function of Mach number shows that some scatter exists. This scatter 
may be attributed to fluctuations in the air supply from the Allis- 
Chalmers compressor. The presence of this problem and methods under- 
taken to alleviate it have been discussed previously in Section 3.1. 

This scatter of data is also present in the curves of as a function 
of Mach number given in Figure 17, but to a lesser extent. 

In order to use the results of the dummy rotor tests in later data 
reduction programs, it was necessary to obtain equations representing 0 
and in the form 

= Constant [f^ (clearance) ][ f^ (Mach number)] 

The function f^ was obtained by first taking the ratio of o? ^ to 
maximum at each Mach number, for each clearance. The values of this 
ratio were then plotted versus clearance. From this plot an average curve 
was obtained, and its equation was computed by the least squares technique. 
The resulting equation expressed °<^/ °^ max as a dimensionless function 
of clearance. 

The function f^ was obtained by computing the ratio of to the 

value of cK ^ at a Mach number of 0.5, and plotting this ratio versus 
Mach number for each clearance. Again, an average curve was determined, 
and its equation was computed. This equation, again dimensionless, ex- 
presses CX^/ 5 ) as a function of Mach number. 

The appropriate scaling factor was determined for each test point. 

An average of all of these factors, which did not differ greatly in value, 
was taken as the constant in the above expression. The equation for 0 
was determined in the same manner. 
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The final equations resulting from this procedure were 



°<i= (8|.762)[.982£7 +■• 2695-8 (Cine) " • 65030 (Cine) 2, ] 

[. 9962^ f . 0Z223(>1i) 02990 (H)*] 

<^ = (• 92^73) [i.Oi433- 1.Z*+601 (Clnc)4 ZH.O(,$39 (One) 2 * „ v 

(14) 

- ISO. 3625-8 ( Cl^c) 3 ] 

[•93iH9f • 18292 09202 (/Mj) 2 ] 



Families of curves produced by these equations are given in Figures 18 
and 19, for (ft and respectively* These curves agree well with those 
of Figures 16 and 17 when consideration is given to the approximations 
introduced in the derivations of Equations (13) and (14). For the veloc- 
ity coefficient, ^ , the maximum difference which exists at any point 
between the curves of Figure 16 and those of Figure 18 is 0.0052. A 
similar comparison for o ( ^ yields a maximum difference of 0*06 degrees. 
Because of the manner in which Equations (13) and (14) were derived, the 
most significant discrepancies occur near the extremes of the Mach number 
range . 
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SECTION 4 



ANALYSIS OF OVERALL TURBINE PERFORMANCE 

The turbine overall performance tests, when combined with the results 
of the dummy rotor tests, provided data from which detailed performance 
parameters of the turbine could be evaluated. The specific parameters 
evaluated were: the total-to-s tatic overall efficiency; velocity and 

loss coefficients for the rotor and guide vanes; rotor inlet and outlet 
flow angles; referred speed, moments, and flow rates; the head coefficient; 
and the degree of reaction. Calculations were made assuming no bearing 
losses and also for the case of "maximum" bearing losses. Bearing loss 
values were determined from equations given by Vavra (1), which were 
established from coast-down tests. The analysis was made using a '‘mean 11 
streamline approach, which is defined as utilizing a representative 
streamline along which the measured flow conditions are assumed to exist. 
This streamline was taken as that at the outer discharge radius of the 
turbine for this analysis, and thus all losses are referred to that outer 
radius. 

4. 1 Instrumentation 

The '•active" rotor shown in Figure 7 was installed for this series 
of tests, and the rotor and wooden casing were placed in an aluminum 
test stand. Figure 1 shows the test installation. Alignment of the scroll 
with respect to the rotor was accomplished by the use of the dummy align- 
ment shaft and dial micrometers shown in Figure 20. Precise alignment 
was found to be quite difficult because the massive wooden casing and 
scroll assembly had to be positioned relative to the rotor. Positioning 
was accomplished by means of four adjusting bolts located at the bottom 
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of the casing and by four angle-irons attached to the casing and resting 
on the bearing supports* In spite of the difficulty involved* alignment 
of the concentric aluminum rings which serve as supports for the shrouds 
was obtained to within + 0*001 inch in the radial direction from the 
dummy alignment shaft* These rings were found to deviate from a plane 
perpendicular to the rotor axis by less than 0*002 inch at a radius of 
4.7 inches* 

The power developed by the turbine was absorbed by a Vortec Products 
Company air dynamometer* model 20*075* shown in Figure 1* Lubrication of 
the dynamometer was accomplished by an oil mist generated by air pressure. 
The torque output was measured with a torque capsule having a maximum 
rated capacity of 400 inch-pounds* This capsule utilized strain gages 
arranged on an internal flexure and was additionally instrumented with an 
iron-constantan thermocouple to provide a means of monitoring the capsule 
temperature* The signal generated by the strain gages was read on a Day- 
tronic model 700 digital indicator located in the control room (Figure 10). 

The dynamometer shroud was modified from the configuration used by 
Riley (3) so as to discharge the hot exhaust air radially outward, rather 
than back onto the torque capsule and housing* in order to minimize ther- 
mal effects on the torque capsule* The circular aluminum ring visible in 
Figure 1 was used to accomplish this modification* Additionally, a wooden 
baffle, not shown on Figure 1* was used to protect the capsule from the 
cold turbine discharge air* and to prevent any interference between the 
dynamometer discharge and the right rotor discharge* 

Figure 21 shows the inner face of the torque capsule, the steel quill 
shaft used to connect the dynamometer to the turbine* and the magnetic pick- 
up and six-lobe flux cutter used to measure turbine speed* The turbine 
speed was read from a Hewlett-Packard electronic counter located in the 
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control room (Figure 22) . Speed was regulated by means of a remotely 
controlled motor on the dynamometer itself, which changes the load capac- 
ity of the dynamometer* The remote control unit is shown in Figure 10, 
and the drive motor is shown in Figure 1. 

Lubrication of the turbine bearings was provided by a dry-sump system 
utilizing a gear-driven oil pump and an upper and a lower oil mist injector 
for each bearing. The lubricating fluid used was 1010 jet engine oil. 

The cart containing the pump, oil reservoir, heat exchanger, and system 
pressure gages is shown in Figure 23, and the lubrication system is shown 
schematically in Figure 24. For these tests the previously existing sys- 
tem was changed to that shown in Figure 24 to permit measurement of the 
oil temperature into the bearings and at the point of discharge from both 
bearings. This was accomplished with iron-constantan thermocouples in- 
stalled as shown on the schematic. 

Figure 24 also shows the iron-constantan thermocouples used to mea- 
sure the bearing outer race temperatures. Vibration signals from the 
bearings were detected by a piezo-quartz accelerometer mounted on the up- 
per left bearing cap and monitored on a Panoramic vibration analyzer manu- 
factured by the Singer Metrics Company. The accelerometer is visible in 
Figure 8, and the analyzer is shown in Figure 22. 

The pressures and flow rates were measured in the same manner as for 
the dummy rotor test previously described. The signals from the thermo- 
couples were initially read from the Brown potentiometer shown in Figure 
13. During the latter part of the testing period, however, they were 
read from a new automatic data acquisition system manufactured by B & F 
Instruments, Inc. This system is shown in Figure 25. 
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4.2 Test Procedures 

Test runs were conducted at five different axial clearances: 0.081, 
0.061, 0.041, 0.024, and 0.015 inches. For each clearance, runs were 
made at total inlet to static discharge pressure ratios of 1.2, 1.3, 1.4, 
1.5, 1.6, and 1.7. At each pressure ratio, turbine speed was varied from 
the minimum speed (maximum load condition) for that pressure, to the 
highest speed attainable, in increments of about 1,500 rpm. Maximum 
speed was restricted to 18,000 rpm for reasons of safety. 

Measurement of the torque at zero speed was also made for each pres- 
sure ratio at each clearance. This was accomplished by engaging the dyna- 
mometer rotor lock and measuring the resulting torque. 

As turbine speed was increased during the first run, a point was 
reached at which the bearings overheated and started to seize, causing 
a rapid drop in speed. The turbine was dismantled and carefully examined. 
No apparent damage had been incurred by the ball bearings. The oil sys- 
tem was flushed and refilled, and the unit was reassembled. This condi- 
tion repeated itself in several subsequent runs, occurring whenever the 
bearing temperature reached approximately 135°F. Additionally, consider- 
able amounts of oil were being lost from around the shaft and bearing 
caps at higher speeds. Apparently neither of these conditions had ex- 
isted in previous tests. 

Since the bearings themselves did not appear to be at fault, ex- 
tensive tests were made to determine whether changing the amount of lubri- 
cating oil flow would alleviate these problems. It was found that vary- 
ing the oil system pressure from 10 pounds per square inch to 50 pounds 
per square inch would reduce the bearing temperature a maximum of six 
degrees, but this variation did not significantly reduce leakage. A 
solution was ultimately found to both problems by carefully positioning 
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the oil injectors so that they discharged directly into the balls of the 
bearing and did not strike the race and by slightly enlarging the injector 
orifices so as to increase the mass flow rate of the oil. The bearing 
temperatures could then be kept below a maximum of about 125 °F at 18,000 
rpm, and oil leakage was negligible. 

The fact that the bearings experienced a maximum temperature in- 
crease of about 50° over the range of speeds tested, whereas in previous 
tests the increase had been only about half as much (7), created a ques- 
tion as to the applicability of the bearing loss information obtained 
from the coast-down tests conducted by Vavra (1). It was therefore decided 
to attempt to evaluate the bearing losses by measuring the oil flow rate 
and the amount of heat absorbed by the oil as it flowed across the bear- 
ings. If the constant pressure specific heat, C^, of the oil were known, 
the energy lost in the bearings could be obtained from the relation 



BFM « 



Won Cp oil bCAClO^S 



XjJ 



(15) 



where BFM is the bearing friction moment in foot-pounds and xo* is the 
turbine rotational speed. 

Numerous tests were made in order to calibrate the mass flow rate of 
the oil with the system pressure, and the lubrication system was modified 
to include the previously discussed thermocouples ahead of and after the 
bearings. Flow rate calibration was accomplished by weighing the amount 
of oil discharged by the injectors on each side in a given time interval. 

In attempting to evaluate bearing losses by this technique, the re- 
sults obtained from several test points indicated losses on the order of 
20 per cent of the minimum losses predicted by Vavra (1). In fact, it 
was discovered that the exposure of the inlet lines, the bearing block it- 
self, and the sumps to the cold discharge air sometimes produced cooler 
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temperatures on the discharge side of the bearings than existed ahead 
of the injectors » Since insulation of the oil lines and the bearing 
blocks to the extent necessary to overcome this deficit was not feasible, 
this attempt to determine the bearing losses was abandoned. The oil 
temperatures were recorded as data for future use 9 however, and appear 
in Appendix R, 

The dynamometer torque capsule was calibrated both before and after 
each day 8 s runs in order to permit consideration of any variations during 
the course of the runs. Calibration was accomplished by use of the de- 
vice shown in Figure 26, An aluminum plug was fabricated to fit into the 
torque ring of the capsule. Torque was applied by loading the weight 
pans on the moment arm, which was attached to the plug, Riley (3) at- 
tached the moment am directly to the inlet side of the dynamometer; 
thus, a moment was introduced due to the axial distance between the 
capsule and the point of application of the load. The apparatus used in 
these tests reduced this potential source of error almost entirely. 

Loads of 0 to 400 inch-pounds were applied to the capsule in 100 
inch-pound increments. The span of the digital readout was set at 0 to 
1000 counts, the latter approximately corresponding to 400 inch-pounds • 
Calibration readings were taken for both increasing and decreasing loads 
in order to evaluate the amount of hysteresis present in the readout. 

As each load was applied, the capsule housing was tapped gently with a 
mallet to overcome any drag in the capsule itself and to eliminate errors 
caused by lag in the servomechanism of the digital counter. This pro- 
cedure is valid, because a like amount of vibration was observed to exist 
under the dynamic load conditions of an actual test run. Using this 
technique, the amount of hysteresis present between the loading and un- 
loading cycles was found to be negligible. 
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The torque readouts for calibrations performed both before and after 
the test runs were found to vary linearly with applied torque 0 However, 
the data obtained after the run were generally higher, by 2 to 3 counts 
at zero load, to 13 to 15 counts at 400 inch-pounds, than the original 

y 

calibration,. This was felt to be due to the thermal effects of the higher 
capsule temperatures which existed during the run* A thermocouple was 
installed in the capsule in an attempt to establish some correlation be- 
tween this temperature rise and the torque count increase* Tests in which 
the capsule was heated from the outside with a strip heater were incon- 
clusive, because similar increase in torque counts for given torque load- 
ings could not be produced, even though the temperatures indicated by the 
thermocouple were considerably higher than those recorded during the actual 
test runs. 

When the turbine is in operation, the capsule is heated from within 
by heat from the dynamometer bearings * This is modified by the dynamometer 
water cooling system and by the cold turbine discharge air* A method of 
reproducing the interaction of all of these conditions in static calibra- 
tion tests was not found, and so the torque calibration data ultimately 
used represented an average of the data obtained before and after the test 
runs. This accounts for some increase due to thermal effects and is felt 
to be a reasonable representation. In some instances, the capsule tempera- 
ture was observed to increase immediately following shut-down, while at 
other times the reverse effect occurred. It is felt that more representa- 
tive temperature data could be obtained by relocating the thermocouple on 
the flexure closer to the strain gages. 

The compressor which provides house air failed after two clearances, 
0.081 and 0.061 inches, had been tested. This compressor provides air to 
the speed controls of the Allis-Chalmers compressor and also provides the 
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air which is used as a reference pressure for the water manometer board. 
The remaining three clearances were tested using a smaller auxiliary 
compressor, which introduced considerable fluctuations into the house air 
system. 

4.3 Theoretical Analysis 

The known quantities necessary for the turbine analysis as performed 
in this project are* mass flow rate; total inlet temperature and pres- 
sure; static inlet pressure; turbine speed and torque output; atmospheric 
pressure and temperature; and the absolute rotor inlet flow angle and 
velocity coefficient based on the results of the scroll and guide vane 
analysis • 

The working fluid is assumed to undergo an adiabatic process, since 
the quantity of heat passing through the casing is negligible compared to 
the enthalpy change undergone by the fluid as it passes through the tur- 
bine. It is additionally assumed that the fluid behaves as a perfect 
gas. 

The flow is assumed to have axisymmetric stream surfaces; thus, in 
effect, the analysis is one-dimensional, based on the determination of 
the flow properties along a particular streamline lying on a chosen sur- 
face of revolution. Vavra (8) states that the major portion of the rotor 
losses is due to secondary flow phenomena, tip leakage losses, and the 
loss that occurs because of the motion of the blade tips in the wall 
boundary layers. Thus is would appear that the majority of the losses 
occur near the shrouds, and the representative streamline was accordingly 
chosen as that existing at the outer radius of the rotor discharge. It 
was felt that this approach would yield a better correlation of results 
than that obtained from use of a calculated mass-averaged mean radius at 
the rotor discharge (7). 



43 



It is assumed that the relative flow angle p calculated from the 



blade dimensions, is very nearly equal to the actual relative flow angle 
at the outer discharge radius (where ^ assumes its maximum value) . This 
implies that the so-called slip factor may be neglected., Verification of 
the validity of this assumption is given by the flow surveys taken at the 
rotor discharge by Riley (3) . Accordingly, j3 was calculated from 
Equation A2(4) of (1) and found to be equal to -69.85 degrees at the outer 
radius. This value was used in all calculations. 

The thermodynamic process of a fluid passing through a turbine is 
given in Figure 14. The resulting flow velocities are shown in Figure 15. 
The derivations which follow are based on these two diagrams. 

Referring to Figure 14, for a reversible, adiabatic process, 



an isentropic expansion through the turbine to p^* The isentropic tempera- 
ture drop which would occur for such an expansion is 



Y-i 




(16) 



where T£ f 9 represents the static temperature at the rotor discharge for 




(17) 



and may therefore be expressed as 




(18) 



Similarly, for an isentropic expansion across the inlet guide vanes 




V-i 



(19) 



and 




( 20 ) 
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The isentropic, or theoretical degree of reaction, r, represents 
that portion of the overall isentropic enthalpy drop which occurs in the 
rotor, or 



r = V- T * 



Therefore, 



r = 1 - 



aT: 






( 21 ) 



e varus 



*T Cs 



(22) 



and, from Equations (18) and (20), 

[l- ( -Pi / Pi. ) ^ J 



r - L- 



(23) 



r i - L?z/pu)*r j 

The total-to-statxc efficiency of the turbine is defined as the 
ratio of the work actually produced by the turbine to that work which 
could ideally be accomplished by isentropic expansion from total inlet 
to static outlet conditions o For this case, the efficiency may be ex- 
pressed by the ratio of the corresponding temperature drops, 



4 T w ^ h T is 



(24) 



The change in total enthalpy across the turbine is 



= c, aT„ = < 25 > 

where M is the net torque produced, Wis the rotational speed, and W 
is the mass flow rate* Therefore, may be expressed ass 



aT w 



M w 

"J Cp W 



(26) 



The efficiency given by Equation (24) was computed both for the 
case where bearing friction was considered and for the case where it was 
not. Bearing friction was calculated based on the results of coast-down 
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tests performed by Vavra (1), The friction horsepower is expressed by 



HP f = - 0.6 + (N/iOoo) (.17143) (27) 

for N < 10,500 rpm, and by 

HPf - - 0-6 f (N /l000)CnN5)-(‘i.S98 ^lo' 3 J(^/iOOO - 1 0.5) 2, ( 28 ) 

for N > 10,500 rpm. The bearing friction moment is then given by 

m. {M) 

* XOr 

The net moment is then the sum of T, the moment measured by the dynamo- 
meter, and the bearing friction moment: 

M - T + Mf (30) 



The velocity C n corresponding to the isentropic enthalpy drop 
i u 

through the turbine is given by the expression 

Co Z 5 2gJc P &T l s (3i) 



and the peripheral speed of the rotor at the rotor inlet is 



U| 



TTNiDi 

720 



(32) 



where is the rotor inlet diameter in inches « By Equations (31) and 
(32), the isentropic head coefficient is given by 



k i5 = (co/u.f <»> 

The absolute velocity is expressed by Equation (4). An alternate 
expression, obtained by substitution of Equation (23) into Equation (4), 
is 
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Vi = <t> J 2g T c f Q-r) 



(34) 



The temperature may then be obtained from 



Tx - T t o - 



vj 



2.g Jc P 



and M^, the Mach number at the rotor inlet, is 



(35) 




V i 

2-gTCpir 



(36) 



The results of the dummy rotor tests can now be used by introducing 
Equation (14) to obtain ^ by an iterative procedure which is discussed 
in Appendix B. Equation (12) gives the nozzle loss coefficient as 



'fw = 1 - 9* *- < 37 > 

Values for o^y can now be taken from Equation (13) to determine the 
peripheral and meridional velocity components of V^. 

Vui - Vi si n (38) 

Vrm - Vi COS (39) 



From Figure 15 it is apparent that the peripheral component of the 
relative rotor inlet velocity is 



W 



Ul 



Vui - Ui 



(40) 



and 



is then found from 






"H Wui 2 



(41) 



47 



The relative inlet flow angle is 

(3i - sin 1 ( Vul j U ) 

At this point all conditions at the entrance to the rotor are known. 



It is assumed that only the kinetic energy V t /2gJC_ is useful in 

ml P 

the rotor (since it has radial blades), and that the incidence loss at 

2 

the rotor inlet is taken as /2gJCp (8) • Thus, the effective rela- 

tive velocity head is reduced, and static conditions at the rotor inlet 
are really those at T^* on Figure 14 * 

The equivalent state point E shown on Figure 14 is established from 
the energy equation for a rotating system in the following manner (6): 



H*. - Octant = Kt + - -JPT ~ ^ f 



Wz _ Rz* 
2-gT 2gT 



(43) 



where H is the relative total enthalpy. This equation may be rearranged 
R 

in the form 



He 




fUi x - u 

V HT 





H J 



(44) 



where the left side is equal to H » the total enthalpy at the state point 
E» Equation (44) clearly shows that state point E is the proper refer- 
ence point from which to view the expansion process through the rotor in 
an analogous manner to that of the stator. The peripheral speed at the 
outer discharge radius, is given by Equation (32) with replaced 

by outer diameter of the rotor discharge. 

The work output of the turbine, represented by the change in total 
enthalpy, may be represented by Euler's turbine equation (6): 



A Hy/ 



Cp aT w 



Vi Vul - Uzc Vuzc 

gT 



(45) 
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where U^q and V^q occur at the outer discharge radius- Manipulating 
this equation and Equations (24) and (31), an expression for the peri- 
pheral component of the absolute discharge velocity, V^q * is 



Vu 



1.0 






Ri 


’ Vui 


- — 1 


' c. y-- 


Rio 


Ui 


2- \ 





(46) 



Using the velocity diagram of Figure 15, it can be shown that the 
relative velocity at the outer discharge radius, ^£0 * ma ^ ex P resse< ^ 



by 



W 20 _ ( Vu *>/ \Ji) - (R a ,/ R x ) 

Ux Sin ^ 2.0 

where, in keeping with the stated assumptions. 



(47) 



- - 69. 85 ° 



(48) 



If the divergence of the pressure lines on the entropy diagram of 
Figure 14, in the direction of increasing entropy, is neglected, there is 

Ti* - TV = T/ - V (w) 



By Equations (21) and (31), the right side of the expression is equiva- 
lent to 



Ti'- Jz - r 



(50) 



2-gTCp 

From Figure 14, the left side may be expressed in terms of the velocities 



V-Tz' 



l^2ot6 



+ 






Vr 



mi 



(51) 



^8 ^" C P 2$"Tcp 

Manipulation of Equations (50), (51), and (32) yields a relation for 



W 20th’ the theoretical relative velocity at the rotor discharge; 



W: 



2 Oih 



Vl 



W - 1 



(52) 
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Using Equations (47) and (52) it is now possible to define a velocity 
coefficient Y f° r the rotor, which is analogous to the velocity co- 
efficient (f> for the scroll and guide vanes given by Equation (2), 




(53) 



The rotor loss coefficient, which is a measure of the kinetic energy loss 
through the rotor, is, then 



'f* = 1 - f * 



(54) 



The meridional velocity component of the absolute discharge velocity can 
be expressed by 



Vnuo - W20 c -°5 ^2.0 



(55) 



from Figure 15. The absolute discharge flow angle is then 



20 ~ t&n Vu zo I 



(56) 



If consideration is given to the case where bearing losses are as- 
sumed to exist, the resulting efficiency will obviously be different 
from that for no bearing losses. Examination of the foregoing analysis 
shows that, correspondingly, different values for V m20* W 20* ^ 

and (j? ^ will result. 

Thus, from the measured items of test data specified at the begin- 
ning of this section, the following performance parameters are obtained 
by the analysis technique: flow rate, head coefficient, degree of re- 

action, absolute and relative rotor inlet flow angles, velocity and loss 
coefficients for the scroll and guide vanes, the net torque produced, the 
absolute discharge flow angle, the rotor velocity and loss coefficients, 
and the total- to-static efficiency. 
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In this experiment the above parameters are referred to standard 
air in accordance with the NASA method, as described in (1). This method 
utilizes a total inlet pressure of 14.7 pounds per square inch absolute, 
a total inlet temperature of 518.7 degrees Rankine, a specific heat ratio 
of 1.4, and a constant pressure specific heat of 0.24 BTU per pound mass- 
degree Fahrenheit. The referral parameters are 



and 



5 5 Pto/ 1*1.7 
0 S tv/ SI 8. 7 



£ 



0 . 8‘o 




Vi I 

2(Y-l) 



(57) 

(58) 

(59) 



The factor £ , which corrects the flow rate through the turbine for 
varying values of Y , is that given by Vavra (1) and is supposed to give 
a better correlation than the expression for € used by NASA. 

The referred speed is then given by 



N 



ref. 




(60) 



The referred flow rate is 



Y 7 ref, r V/ -f- ) < 61 ) 

Finally,, the referred moment is 

M ref E < 62 > 
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4.4 Data Reduction: Programs RADIAL and NOSPD 

The computer program RADIAL, listed in Appendix B, determines the 
turbine performance parameters described in Section 4*3 from the test data. 
In general, the specific techniques used by the computer in these computa- 
tions are apparent from the program itself. However, Appendix B describes 
each subroutine briefly and elaborates on specific computational pro- 
cedures which are not immediately obvious. 

Program NOSPD, listed in Appendix C, was written to compute the de- 
sired performance parameters from the locked rotor test data. 

4.5 Discussion of Results 

The input data and results of program RADIAL are listed in Tables B4 
through B13, with input and output for a given clearance appearing in ad- 
jacent tables. The input data and results of program NOSPD appear in 
Table C2. 

Figure 27 shows the degree of reaction, r, as a function of head co- 
efficient, k^ g , for the clearances tested in this analysis as well as 
those of Riley (3) and Vavra (1). For each clearance, the degree of re- 
action was found to be a unique function of head coefficient, independent 
of turbine pressure ratio. This is in agreement with the findings of 
both Riley and Vavra. Riley further concluded that the degree of re- 
action is independent of axial clearance. The results of the present tests 
verify this for all clearances except for 0.081 inch. At this clearance, 
the values of r obtained were found to be approximately 1 to 1.5 per cent 
lower than those for the other clearances. This is shown on Figure 27, 
where the resulting curves for all clearances except 0.081 inch are rep- 
resented by that for 0.061 inch, and the curve for 0.081 inch is plotted 
separately. This discrepancy implies that, for a given overall pressure 
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ratio, the pressure p, can no longer be maintained at a clearance of 
0*081 inch. 

It may be seen from the representative curve of Figure 27 that r 
decreases with increasing k^ g (decreasing speed), from a value of ap- 
proximately 0.75 at a k 0 of 1.05 to 0.30 at a k. of 6. This drop oc- 

is is 

curs because, with decreasing speed, a smaller pressure drop is required 
to pass the flow through the rotor. 

The degree of reaction curves obtained in earlier tests by Vavra (1) 
and Riley (3) are superimposed on the curve of Figure 27. Vavra states 
that the value of r at design k^ of 2.3 should be approximately 0.486, 
as compared with the value of 0.395 obtained by him and with Riley’s value 
of 0.402. The current tests yielded a value of 0.440 at design k^ g , 
except for a clearance of 0.081 inch, for which the corresponding value 
of r was 0.434. This apparent increase in degree of reaction is due pri- 
marily to a more accurate determination of p^ than was used by either 
Vavra or Riley. Boshoven (4) relocated the shroud pressure taps, which 
measure the pressure p^, to a position axially outward in the annulus 
around the periphery of the shroud. This was done to reduce the effects 
of local perturbations on the pressure readings and greatly improved the 
uniformity of the pressure distribution around the shrouds. This im- 
provement in the accuracy of p^ yields a more accurate value of r since 
degree of reaction is directly dependent on p^. 

The degree of reaction remains lower than the design value because 
the guide vane exit area is too small to permit discharge of the design 
mass flow (1). 

Significant plots resulting from the overall turbine efficiency cal- 
culations are given in Figures 28 through 31. Figure 28 depicts turbine 
efficiency as a function of referred speed for the case of no bearing 
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losses and a clearance of 0.081 inch. Figure 29 is a corresponding plot 
of efficiency versus head coefficient. Figures 30 and 31 represent these 
same respective parameters, also at a clearance of 0.081 inch, for the 
case where bearing losses are considered. 

Efficiency was found to be essentially independent of clearance over 
the range of clearances tested. Figures 28 through 31 are therefore rep- 
resentative of the results for all clearances between 0.015 and 0.081 inch. 

Values of turbine efficiency were obtained from measurements of the 
torque at the turbine shaft made with the air dynamometer and from the 
bearing loss values as determined from the results of coast-down tests 
conducted by Vavra (1). Because of the previously discussed problem of 
excessive bearing temperatures and seizing experienced during these tests, 
the bearing losses are probably greater than the maximum losses deter- 
mined by Vavra. However, since time did not permit a repetition of the 
coast-down tests, Vavra' s equations for maximum bearing losses [Equa- 
tions (27) and (28)] were used. 

A comparison of Figures 28 and 29 and of Figures 30 and 31 shows 
the advantage derived from representing efficiency as a function of head 
coefficient instead of referred speed. Use of the head coefficient ef- 
fectively reduces a family of curves to a single curve. However, if the 
bearing losses are ignored, as in Figure 29, the resulting plot depends 
significantly on pressure ratio. In this figure only those curves rep- 
resenting the pressure ratios for maximum and minimum efficiencies are 
shown to avoid clutter. 

Figure 31 shows that the efficiencies are nearly independent of pres- 
sure ratio in the case where allowance is made for bearing losses. If 
these bearing losses were accurately known, it would appear that the in- 
fluence of pressure ratio would be completely eliminated from a plot of 
efficiency versus head coefficient. 



54 



In general, if bearing losses are not considered, a pressure ratio 
of 1.7 yields a maximum efficiency of 81.1 per cent at a of 2.35 and 
an axial clearance of 0.015 inch. Minimum efficiencies are obtained for 
a pressure ratio of 1.2 and are in most cases 10 to 12 per cent lower 
than maximum efficiencies in the vicinity of the design k Jg of 2.3. 

Vavra (1) reported a maximum efficiency of 81 per cent for the case of no 
bearing losses at a pressure ratio of 1.616, k^ g of 2.37, and clearance 
of 0.062 inch. Riley (3) obtained a corresponding maximum efficiency of 
82.52 per cent at a of 2.26, pressure ratio of 1.678, and clearance 
of 0.027 inch. 

For the case where bearing losses are taken into account, the maxi- 
mum efficiency, ^ , was found to be 84.6 per cent at a pressure ratio 
of 1.7, head coefficient of 2.28, and clearance of 0.081 inch. This 
compares favorably with the design efficiency of 85.5 per cent at a pres- 
sure ratio of 1.7 and head coefficient of 2.3, as given by Vavra (1). 
Riley (3) obtained a maximum efficiency of 86.04 per cent at a pressure 
ratio of 1.678, head coefficient of 2.259, and clearance of 0.027 inch, 
based on the maximum values of bearing losses. 

At minimum speed (maximum k^ g ) for a given pressure ratio and clear- 
ance, the efficiency in general is about 65 per cent. As speed is in- 
creased (k^ g decreased), the efficiency increases to a maximum of about 

84 per cent in the vicinity of a k. of 2.3 and then falls off rapidly 

r s 

to 30 to 40 per cent at a k 0 of about 1.05. Initially, then, the ef- 

1 s 

ficiency increases with speed because, by Equation (26), an increase in 
Ij and decrease in W offset the decrease in torque, resulting in in- 
creased turbine work and efficiency. As speed is further increased, the 
decreasing torque begins to predominate, causing the efficiency curve to 
peak and then rapidly decrease. 
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It was expected that, as the clearance was varied over the range en- 
compassed by these tests, a ‘''critical 81 clearance would be reached beyond 
which turbine efficiency would deteriorate rapidly. Riley (3) observed 
this to occur at a clearance of 0*052 inch. Such an occurrence was not 
observed in these tests over the range of clearances from 0.015 to 0.081 
inch. Thus, it appears that axial clearance has no appreciable effect 
on turbine efficiency over this range of clearances. 

It has been previously stated that Riley (3) was in all likelihood 
operating at a considerably greater clearance than he believed himself 
to be. This is substantiated to a degree by the fact that Vavra (1) did 
not experience this efficiency drop at a clearance of 0.062 inch. Since 
gaps also existed between the rotor and shrouds during Riley f s tests, 
the blades would have left the wall boundary layer at a lesser axial 
clearance than in these tests where the blade and shroud contours were 
matched. Thus, it is reasonable to assume that the critical clearance 
is at some value greater than 0.081 inch. 

That such a drop in efficiency should occur as clearance is increased 
to some critical value is explained in part by Csanady (9). The rotating 
blades M scrape up 88 the shroud boundary layer and produce vortices in the 
vicinity of the tips of the blades which, in turbines, tend to nullify 
the tip vortices themselves. At a particular tip clearance, then, these 
two vortices would neutralize each other, leaving only the passage vortex 
which causes secondary flow losses. This neutralization should produce 
a pronounced increase in efficiency at that clearance at which it occurs. 
For the subject turbine, it would appear that if this condition does 
exist, it occurs at a clearance less than 0.015 inch. 

The efficiency is essentially independent of clearance because, al- 
though the tip vortices become more predominant as clearance increases. 
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there is also a buildup of the boundary layer. This buildup maintains 
a reasonably even relative balance between the tip vortices and the 
"scraped up 55 vortices. As clearance is extended beyond some critical 
value, however, the blades begin to operate near the edge of the wall 
boundary layer, and the countering vortices quickly cease to be effec- 
tive in neutralizing the increasing tip vortices. The result is a pro- 
nounced increase in rotor loss and a corresponding decrease in turbine 
efficiency. It is also worthy of note that, for the test turbine, al- 
though the axial tip clearance is increased, the radial clearance between 
the shrouds and the rotor blades in the vicinity of the discharge remains 
essentially unchanged. 

The rotor velocity coefficient, ^ , is plotted as a function of 

the head coefficient in Figure 32, for the case where bearing losses are 

neglected. The clearance represented is 0.024 inch. Figures 33 through 

37 show the curves of versus k , with bearing losses considered, for 

is 

each of the five clearances tested. Again, it may be seen that if bear- 
ing losses are ignored, there is considerably more variation of the 
curves with pressure ratio than in the case where bearing losses are 
considered . 

The only variation in rotor velocity coefficient with axial clear- 
ance which can be established is that values of ^ for a clearance of 
0.081 inch seem to be generally higher than for all other clearances. 

This variation is evident for both bearing loss cases, with values of 
being generally 10 to 15 per cent higher at a clearance of 0.081 
inch. Because of the negligible effect of axial clearance and the doubt- 
ful significance of the case where bearing losses are neglected. Figure 
32, for a clearance of 0.024 inch, is the only curve included for the 
case where bearing losses are not considered. 
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Figures 33 through 37 show that maximum values of do not always 

occur at the higher pressure ratios , although Yl is more nearly con- 
stant over the range of values for the higher pressure ratios. 

At the design of 2.3, a value for of 0.960 was produced at 

a pressure ratio of 1.2 and clearance of 0.081 inch. The maximum value 

of was found to occur at the same pressure ratio and clearance and 

at a k„ of 3.35. Its value is 0.990. 
is 

It can be seen from the curves for given in Figures 33 through 

37 that remains relatively constant over a range of k^ g between 

about 2.3 and the maximum value. As k, is decreased below 2.3. V, 

IS 9 XL 

first decreases, with a minimum occurring at a k. of 1.7 to 1.9. It 

is 

then increases to a peak at a k^ g of about 1.5, and falls off sharply 
as the head coefficient is further decreased. Examination of the tabu- 
lated performance data in Appendix B shows that the relative rotor inlet 
flow angle |8i is positive at values of k^ greater than about 2.3, so 
that the relative flow strikes the pressure side of the blades. This 
produces the accompanying incidence losses, and further losses occur as 
the flow is turned through a relatively large angle to follow the flow 
channel. 

As U- is increased (k. decreased), /3i becomes negative, so that im- 
pingement of the flow on the blades acts opposite to the direction of ro- 
tation. The incidence losses are not significantly changed; however, the 
flow must now be turned by a considerably smaller amount in order to 
follow the flow channel. Thus, the resultant sum of these three losses 
appears to be less, and an increase in ^ occurs at a k^ of about 1.5. 
Below this point, as becomes increasingly negative, the relative flow 

becomes almost perpendicular to the blades. Thus, the advantage is lost, 
and the rotor efficiency deteriorates rapidly. 
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A plot of referred moment with bearing losses included, as a function 
of referred speed, is shown in Figure 38 for data obtained at a clearance 
of 0.081 inch. The values of the torque at zero speed resulting from the 
locked rotor tests are also shown. The effects of clearance on this para- 
meter were also found to be negligible over the range of clearances from 
0.015 to 0.081 inch, there being only slight random variations present in 
the resulting data. Therefore, Figure 38 is a representative plot for 
all clearances. 

As expected, the maximum moment of about 16.1 foot-pounds was found 
to occur at zero speed and at a pressure ratio of 1.7. For the maximum 
speed of about 18,000 rpm, this pressure ratio yielded 6.3 foot-pounds 
for zero bearing losses and 6.8 foot-pounds for maximum bearing losses. 

The lowest pressure ratio, lo2, consistently yielded the lowest values 
of moment. The moment curves of Vavra (1) and Riley (3) are almost 
identical to those obtained in these tests, except for the results obtained 
by Riley beyond the "“critical™ clearance. 

In general, the torque produced by the flow energy at a given pressure 
ratio is a maximum at zero speed and decreases with increasing turbine 
speed. This may be explained by the law of angular momentum, which, by 
reference to the velocity diagram of Figure 15, can be written as 



also increases with speed. Thus, it is seen that the torque produced, 
M, decreases with increasing turbine speed. 






(63) 



For a given pressure ratio, is nearly constant with speed, and the 

« 

mass flow rate W decreases with increasing speed. From Figure 15, V £ 
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After completion of the test runs previously described, an 
informal trial run was made at a clearance of 0*100 inch in an attempt 
to determine the critical clearance* Although not conclusive, the 
results indicated no significant change in the values of the perform- 
ance parameters* This implies that a reasonable overall efficiency may 
be maintained at clearances which might considerably exceed those tested 
in this analysis* Since these additional tests were not complete, the 
resulting data are not included in this report. 
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SECTION 5 



CONCLUSIONS AND RECOMMENDATIONS 

At a pressure ratio of 1,7, a clearance of 0.081 inch, and a head 
coefficient of 2.28, a maximum total-to-static overall efficiency of 84.6 
per cent was obtained for the case where bearing losses were considered. 
The efficiency, referred moment, and degree of reaction were found to be 
independent of clearance over the range of clearances between 0.015 and 
0.081 inch. A clearance tolerance of this sort is very important in 
small turbines and in turbines where thermal expansion must be considered. 
It is recommended that additional testing be performed in order to deter- 
mine the precise clearance beyond which turbine performance deteriorates. 
Tests might also be conducted in which the clearance at the rotor dis- 
charge is varied in addition to the tip clearance, although this would re- 
quire cutting the shrouds or the manufacture of additional shrouds with 
various annular openings. 

The results of this investigation show that computations based on 
flow conditions at the outer discharge radius are simpler and more 
straightforward than those obtained by referring the losses to a calcu- 
lated mass-average radius, and that this simplicity does not compromise 
the quality of the results. On the contrary, a correlation of the rotor 
velocity coefficient is obtained by this method, whereas previous ef- 
forts to obtain meaningful loss parameters were not successful (7). 

There is strong evidence that the bearing loss data used in this 
investigation are no longer applicable. Correct data should be obtained 
either through new coast-down tests or through some implementation of 
the lubricating oil temperature differential across the bearings, as pre- 
viously described. The latter method would be highly preferable since it 
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would provide accurate loss information on a current basis and would cor- 
rectly reflect any change in the mechanical condition of the bearings. 
Successful implementation of such a method would probably be contingent 
on a complete redesigning of the test rig, so that the oil lines are 
routed outside of the discharge flow as much as possible and are thoroughly 
insulated. 

Since any inaccuracies in torque measurement introduce considerable 
error into the calculations as performed herein, additional effort is re- 
quired in order to determine how to properly compensate the torque capsule 
output for thermal effects. An alternative solution would be to utilize 
a force (reluctance) capsule and a moment arm arrangement. This capsule 
could be remotely located and shielded so as to be free of the effects 
of the dynamometer bearing heat and the cold turbine discharge air. 

It is recommended that the modifications described by Vavra (1) be 
performed in order to improve the component matching of the turbine and 
that further tests be carried out in order to determine the effects of 
these modifications on the performance parameters. Specifically, the 
enlarging of the guide vane discharge area or the reduction of the rotor 
discharge area would ultimately increase the relative rotor discharge 
velocity which would reduce the peripheral velocity V^- Greater 
turbine work, and thus higher efficiencies, should result. 

It is felt that the turbine test rig is now in much better operating 
condition than for any previous tests. The instrumentation has been proven; 
and, for the first time, the axial clearances may be accurately set. 
Accordingly, further testing to verify and expand on the performance para- 
meters obtained in this investigation is recommended. 
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FIGURE 2 

DUMMY ROTOR TEST INSTALLATION 
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FIGURE 3 

WOODEN SCROLL INSERT 
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FIGURE 5 

CROSS-SECTION OF TURBINE 
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FIGURE 6 

DUMMY ROTOR AND FLOW STRAIGHTENERS 
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FIGURE 8 

LEFT TURBINE DISCHARGE 
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FIGURE 9 

TURBINE INLET PIPING, FLOW CONTROL VALVES 
AND FLOW MEASURING ORIFICE 




FIGURE 10 

REMOTE CONTROL FOR FLOW REGULATION, 
VORTEC REMOTE THROTTLE CONTROL, 

AND DAYTRONIC STRAIN GAGE DIGITAL INDICATOR 
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FIGURE 11 

SCHEMATIC OF RADIAL TURBINE INSTRUMENTATION 
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FIGURE 14 

TEMPERATURE-ENTROPY DIAGRAM FOR RADIAL TURBINE 
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FIGURE 15 

VELOCITY DIAGRAMS FOR RADIAL TURBINE 
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FIGURE 17 

ABSOLUTE ROTOR INLET FLOW ANGLE VS. MACH NUMBER 
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FROM DERIVED EQUATION 



80 




FIGURE 19 

ABSOLUTE ROTOR INLET FLOW ANGLE VS. MACH NUMBER 
FROM DERIVED EQUATION 
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FIGURE 21 

FLUX CUTTER, QUILL SHAFT, 

AND DYNAMOMETER TORQUE CAPSULE 




FIGURE 22 

ELECTRONIC SPEED COUNTER 
AND VIBRATION ANALYZER 
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FIGURE 23 

FIVE INCH INLET PIPE AND BEARING LUBRICATION UNIT 
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FIGURE 24 

SCHEMATIC OF TURBINE LUBRICATION SYSTEM 
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FIGURE 26 

TORQUE CAPSULE CALIBRATION APPARATUS 
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FIGURE 27 

DEGREE OF REACTION VS. HEAD COEFFICIENT 
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FIGURE 28 

OVERALL EFFICIENCY (NO BEARING LOSSES) VS. REFERRED SPEED 
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FIGURE 29 

OVERALL EFFICIENCY (NO BEARING LOSSES) VS. HEAD COEFFICIENT 
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FIGURE 30 

OVERALL EFFICIENCY (WITH BEARING LOSSES) VS. REFERRED SPEED 
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FIGURE 31 

OVERALL EFFICIENCY (WITH BEARING LOSSES) VS. HEAD COEFFICIENT 
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FIGURE 32 

ROTOR VELOCITY COEFFICIENT (NO BEARING LOSSES) VS. HEAD COEFFICIENT 
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FIGURE 33 

ROTOR VELOCITY COEFFICIENT (WITH BEARING LOSSES) VS. HEAD COEFFICIENT 

CLEARANCE 0.081 IN. 
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FIGURE 34 

ROTOR VELOCITY COEFFICIENT (WITH BEARING LOSSES) VS. HEAD COEFFICIENT 

CLEARANCE 0.061 IN. 
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FIGURE 35 

ROTOR VELOCITY COEFFICIENT (WITH BEARING LOSSES) VS. HEAD COEFFICIENT 

CLEARANCE 0.041 IN. 
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FIGURE 36 

ROTOR VELOCITY COEFFICIENT (WITH BEARING LOSSES) VS. HEAD COEFFICIENT 

CLEARANCE 0.024 IN. 
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FIGURE 37 

ROTOR VELOCITY COEFFICIENT (WITH BEARING LOSSES) VS. HEAD COEFFICIENT 

CLEARANCE 0.015 IN. 
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REFERRED MOMENT (WITH BEARING LOSSES) VS. REFERRED SPEED 
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APPENDIX A 



PROGRAM SCROLL 

Program SCROLL computes the losses in the scroll and guide vanes and 
the absolute rotor inlet flow angle from torque, mass flow rate, pressure, 
and temperature data obtained in the dummy rotor tests The program is 
essentially that used by Riley (3) and by Boshoven (4), with several 
changes as specified below. It can process any number of runs, with a 
maximum of ten sets of data per run. A listing of program SCROLL is 
given in Table Al. The input and output values which constitute the 
results of the eight sets of dummy rotor test data obtained are pre- 
sented in Tables A2 through A9* 

Al Main Program 

The value of the specific gravity of mercury, G^g* at room tempera- 
ture t , was determined from 
rm 

0-36303 x| o‘ 3 ) tr* (Al ) 

for temperatures between 32°F and 113°F. The factor for converting 

2 

inches Hg to lb/ft is 

Cfi - ( Gto g ^13.54^ (A2) 

The value of the specific gravity of water at room temperature is 

Ghjo s 0. 9993763 4* (l‘04e^76 - (/^53 186 Xfo+) tr* < A3 > 

The specific gravity relations were obtained from the tabulated data of 
Reference 10. . These equations are slightly different from those used 
previously. 
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A2 Subroutine TEMP 



The relations for the evaluation of the temperatures in this sub- 
routine from the thermocouple readings in millivolts are different from 
those previously used because of the substitution of iron-constantan 
thermocouples for the chromel-alumel type. Thus, for iron-constantan 
thermocouples producing a millivolt output V, and using a cold junction 
reference temperature, t c j» the temperature relations are 

t = tc] + o.m + -*>8.77 V - o.q5i8 \j z <A4> 

for t s loo 0 P 

t = tcj + I.25Z 4- V - 0. i $55 V z < A5 > 

for 100 'F <t f 200 a F 

Equations (A4) and (A5) were obtained from data tabulated in Reference 

11 . 



A3 Subroutine FLOW 

The evaluation of the flow rate using the measured data from the 
particular flow orifices employed in the radial turbine test setup was 
carried out in accordance with a method developed originally by Vavra (1), 
after the techniques given by Stearns, et al. (5). The values thus deter- 
mined by Vavra applied to a sharp-edged orifice of 2.800 inches diameter. 
Before this series of tests was begun, the flow orifice diameter was 
measured and found to be 2.7965 inches. Additionally, the conversion 
factors and the value of the specific gravity of mercury used herein are 
slightly different than those previously used, the latter being 13.54 vice 
13.59. The equations given below are thus those of Vavra (1), with these 
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modifications introduced. Computations were made for the vena contracta 
taps only, although flange tap data were available, because these data 
give a more accurate flow rate (1). 

The relation for the flow rate is 



where : 
C 



W = c X Fr / ' 



(A6) 



- factor dependent on orifice diameter, type of pressure 
taps, and dimensional units 

- area multiplier to account for thermal expansion of 
orifice 



- factor which accounts for compressibility effects 

- Reynolds number correction factor 

- absolute pressure at upstream tap 

- pressure differential across orifice 

- temperature ahead of orifice 
For vena contracta taps, C - 0.9057. 

For a stainless steel orifice: 



1VC 
Ah 



1VC 



oC 


= to + (1.85 *10' s ) (TV-53Z.) 


(A7) 


Yi 


= 1.0 - 0.351 - - 

n v c 


(A8) 


Fr 


= 1.0 + 

X 


(A9) 


X 


= 0.8151 


(A10) 


z 


- 1.9+ (2.M * iO' 3 )(T^- 659.69) 


(All) 



where Equation (All) is valid for air between 50°F and 300°F. 
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The resulting flow rate equation used herein is then given by 



W = 2 - 28 SH « Yl Fr / 



Pi. vc A hive 

Th 



where: W is in lb /sec 

m 



^lYc® ^ ^IVC in® Hg . 

T, is in °R 
4 

and where and are as given previously. 



(A12) 



A4 Subroutine PRESS 

It has been previously mentioned that a regulated reference pressure 
was used on the 96~inch water manometer board. This pressure must there- 
fore be accounted for in computing the absolute static pressure at 

the rotor inlet. This is done in the equation 



~p _ (Tklrvy 4 Hrgfj ) (Grhff) f ( Hjl) ((XHio) 
ri " 13 . 54 



(A13) 



where the quantity (^ atm represents the total reference pressure 

applied to the manometer reservoir. 



A5 Method of Assembly of Input Data 

The items of input data necessary in performing computations with 
program SCROLL are given below in the order in which they must appear 
on the data cards. Specific formats for the numerical entries are given 
in the program listing. 



Card No . FORTRAN 

1 NRUNS 

2 NPTS 

3 DPVC 



Description 

Number of runs (clearances) to be processed 

Number of points (pressure ratios) in run 

Measured pressure differential across flow 
orifice j, from vena contracta taps (in. Hg) 
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Card No » 


FORTRAN 


Description 




PUVC 


Measured pressure upstream of flow measuring 
orifice , from vena contracta taps (in. Hg 
gage) 




P5P 


Measured static pressure at turbine inlet 
(in. Hg gage) 




PATM 


Barometric pressure (in. Hg) 




HREF1 


Measured reference pressure applied to H^O 
manometer board (in. Hg gage) 




HI 


Average static pressure ahead of rotor 
(in. H 2 0 gage) 




SR 


Torque scale reading (lb.) 




TRM 


Control room temperature (°F) 




V4 


Thermocouple reading ahead of orifice (mv.) 




V5 


Thermocouple reading at turbine inlet (mv.) 


4 


TCJ 


Cold junction temperature (°F) 




STARE 


Torque scale tare (lb.) 




B1 


Perpendicular distance between shroud inner 
extremities (in.) 


Cards 


3 and 4 comprise the input data for one test point and are 


repeated for a total 


of NPTS test points for each run. Each new run 


begins with 


card 2. 
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APPENDIX B 



PROGRAM RADIAL 

Program RADIAL computes the referred flow rate, head coefficient, 
degree of reaction, relative and absolute rotor inlet flow angles, velocity 
and loss coefficients for the scroll and guide vanes, referred moment, 
absolute discharge flow angle, velocity and loss coefficients for the rotor, 
and efficiencies from the torque, speed, mass flow rate, and the pressure 
and temperature data obtained during the active rotor tests. The calcula- 
tion technique generally follows that given in the theory of Section 4.3. 

It can process any number of clearances (sets) and any number of pressure 
ratios (runs) per clearance with a maximum of ten different speed settings 
(points) per run. A listing of program RADIAL is given in Table Bl. 

In program RADIAL, use was made of temporary memory disc storage of 
the input data in order to avoid the awkward procedure of subscripting the 
input variables. This procedure permits the processing and printing of 
the input data on a point-by-point basis and the presentation of the results 
in a compact, tabular form. Appropriate comments appear in the main pro- 
gram at the points where this procedure was used. The additional control 
cards required to implement this technique and to obtain necessary addi- 
tional memory core storage space are shown in Table B2. As the table indi- 
cates, these cards are inserted just ahead of the data cards in the as- 
sembled deck. 

The torque calibration data used for all runs are given in Table B3. 

The input and output values resulting from the five sets of test data are 
presented in Tables B4 through B13. 
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B1 Subroutines TEMP and FLOW 



These subroutines perform the same functions as in program SCROLL, 
namely those of obtaining temperatures from the thermocouple readings and 
calculating mass flow rate, respectively « Subroutine TEMP has been ex- 
panded to include calculation of the temperatures of the bearing lubricat 
ing oil at the inlet and at the discharge of each bearing. However for 
reasons discussed in Section 4, these values were not used in this experi 
ment. 



B2 Subroutine PRESS 

Subroutine PRESS establishes the total- to-static pressure ratios 
across the scroll and guide vanes (P^/p^) and across the turbine (P / 
p^) and the ratio P^/P 2 ° An iteration process is necessary in order to 
determine an average value of the total pressure at the turbine inlet, 

P . Three relations are used to accomplish this, these being the gas 
law. 



?. = 

the continuity equation, 

V. = 





w 

f* As 



and the energy equation. 



To 



Tto 



Vo 2 

T c P 



(Bl) 



(B2) 



(B3) 



where is the area of the five-inch inlet pipe and C p is the specific 
heat of air at the temperature T Using T^q as a first approximation 
of an initial value of is found, and from this, by Equation (B3) , 

a new value of Tq. The iteration continues until a difference of 0.01 
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Since 



degrees or less exists between any two successive values of T^. 
velocities are low, the total pressure is then obtained from 

Pt. = + /° v ? - «*> 

Cfi 

B3 Subroutine EDC 

This subroutine computes the mean values of Y and Cp based on the 
arithmetic mean of total turbine inlet and static isentropic discharge 
temperatures* From Reference 3, the variations of Y and with tempera- 
ture are 



Y= Moig - (2* io" ff ) t < B5 > 

Cp = 0.739H3 f ( 3 . >4 *io“*) t 4- (2 * I0~*) t 1 ( B6 > 



These relations were originally obtained from data tabulated in Refer- 
ence 12 and are valid from approximately 40°F to 170°F. 

The values of V and v correspond to the average tempera- 

avg P(avg) 

ture through the turbine based on an isentropic expansion. This average 
temperature, t, is given by 



t - (Tt. - ^5 9. 6?) 



at * cs 
% 



(B7) 



where 



A — Tto 



L- 



(B8) 



Using an iteration process ,, the first approximation of 4 T. is based on 
a value of Y corresponding to T^. The value of t resulting from Equa- 
tion (B6) is compared with the previous value of t in each iteration 
until the difference between any two successive values of t is less than 
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1.0°. Using the final value of t, values of V avg , C p^ avg )» and ^T ig 
are computed from Equations (B5) , (B6) , and (B8), respectively. The 
reference parameter 9, given by 

e= (y.vj/1.4 ) (Ts/sis.t) (B9> 

is also computed in this subroutine. 

The subscript indicating the average values of X and C p will be 
dropped in subsequent discussions. The average values of these quanti- 
ties are assumed in the remainder of the program unless otherwise stated. 

B4 Subroutines TORQ and DYNA 

These subroutines compute the net torque produced by the turbine by 
evaluating the bearing friction moment and the torque measured by the 
dynamometer. 

In subroutine TORQ, the bearing friction moment is evaluated from 
turbine rpm by Equations (27), (28), and (29). The net moment is calcu- 
lated by Equation (30), where the dynamometer torque output is found in 
subroutine DYNA, called from subroutine TORQ. 

Subroutine DYNA determines the torque from the dynamometer calibra- 
tion data obtained for each run. These data are tabulated in Table B3 
for all runs. 

The values of the torque calibration data constitute a one-dimensional 
array with five elements, representing the torque scale reading for loads 
of zero to 400 inch-pounds, respectively, taken at 100 inch-pound intervals. 
Thus, TCD(l) represents the torque scale readout value for the no load 
condition. Using a DO loop with index I, the torque indicator reading 
for the test point being processed, TQ, is compared with the values of 
TCD(I) to determine the interval (I, 1+1) within which TQ lies. Since 
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the calibration curve data are very nearly linear, the torque is computed 
assuming a straight line approximation between TCD(I) and TCD(I-l), by 
the equation 



T - 100(1-1) + 100 



TQ - TCP(r) 
TCD (r + D- tcd (i) 



This procedure was adapted from that used by Riley (3). 



(BIO) 



B5 Subroutine TURB 

The turbine performance parameters described in Section 4 are cal- 
culated by subroutine TURB in a manner which follows the theory previously 
given. No iterations are involved, and the equations used are self- 
explanatory. Therefore, the steps will not be discussed here. The veloc- 
ity coefficient, (p , and the absolute rotor inlet flow angle, are 

provided from subroutine PHIAL, for which TURB is the calling subroutine. 



B6 Subroutine PHIAL 

This subroutine calculates the values of p , o( V^, and ^ 
from the results of program SCROLL. The values of Y and C p used are 
based on the total inlet temperature, rather than the average tempera- 
ture as computed in subroutine EDC, because the values of and o( ^ 
obtained in program SCROLL were computed in this manner. 

The values of (j> and are obtained from Equations (13) and (14), 
for which the entering arguments are CLNC, the axial clearance, and ACH1, 
or , the Mach number at the rotor inlet. Since and T^ are not known 

at this point in the computations, they must be obtained by iteration. 

An initial value of is obtained from Equation (4) by assuming <fi = 1.0, 
and initial values of and are obtained from Equations (35) and (36), 
respectively. Using this value of M^, called ACHIA, and CLNC, a value of 
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0, called PHXA, is obtained from Equation (13). The assumed value of 0 
is reduced in increments of 0.001 until succeeding values of 0 and 0 

Jri 

differ by no more than 0.001. Thus, values of M^, and 0 are estab- 
lished, and c<^ is then computed from Equation (14). The loss coeffi- 
cient for the scroll and guide vanes, is then obtained from Equation 

( 12 ). 

B7 Subroutine REFER 

This subroutine computes the referred speed, referred flow rate, 
and referred moment in accordance with Equations (60), (61), and (62), 
respectively. 

B8 Method of Assembly of Input Data 

In order to facilitate future use of program RADIAL, the order of 



assembly of the input data cards is given. The actual formats for the 
numerical entries may be obtained directly from the program listing. 

All items of data are given in the order in which they must appear on the 
data cards. For convenience, the symbols are redefined. 

Card No. FORTRAN Description 


i 


NSETS 


Number of sets of data (clearances to be 
processed 


2 


CLNC 


Clearance for set 


3 


NRUNS 


Number of runs (pressure ratios) in data set 


4 


TCD 


Torque calibration data for run 
(5 values) . 


5 


NPTS 


Number of data points (speeds) in run 


6 


RPM 


Measured turbine speed (rpm) 




TQ 


Torque indicator reading (counts) 




PUVC 


Measured pressure upstream of flow measuring 
orifice, from vena contracta taps (in. Hg gage) 
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Card No , FORTRAN Description 



P5P 


Measured static pressure at turbine inlet 
(in. Hg gage) 


HI 


Average static pressure ahead of rotor 
(in H^O gage) 


DPVC 


Measured pressure differential across flow 
orifice, from vena contracta taps (in. Hg) 


PATM 


Barometric pressure (in. Hg) 


HREF1 


Measured reference pressure applied to water 
manometer board (in. Hg gage) 


V4 


Thermocouple reading ahead of orifice (mv.) 


V5 


Thermocouple reading at turbine inlet (mv.) 


7 V20 


Thermocouple reading at lube oil inlet (mv.) 


V21 


Thermocouple reading at lube oil discharge 
from right bearing (mv,) 


V22 


Thermocouple reading at lub oil discharge 
from left bearing (mv.) 


TRM 


Control room temperature (°F) 


TCJ 


Cold junction temperature (°F) 



Cards 6 and 7 comprise the input data for one test point and are 
repeated for a total of NPTS test points for each run. Each new run 
begins with card 4, and each new set begins with card 2. 
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END 



/* 

//G0.FTO6FOO1 

//G0.FTO8FOO1 

// 

//G0.SYSIN DD* 

1 i 


DD SYSOUT=A, SPACE=(CYL, (15,1)) 

DD UNIT=SYSDA,SPACE=(TRK,(2,1)) , * 

DCB= (RECFM=V , LRECL=500 ,BLKSIZE=504) 


1 7 


16 COLUMN 72 



TABLE B2 

ADDITIONAL CONTROL CARDS REQUIRED FOR PROGRAM RADIAL 
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TABLE B3 



TORQUE CALIBRATION DATA 



SET 


CLEARANCE 


RUNS 




APPLIED TORQUE (IN-LB) 




(IN) 




0 


100 


200 


300 


400 






1-4 


5. 


261. 


518. 


768. 


1013. 


1 


0.081 


5-6, 
NOS PD 1 


2. 


255. 


510. 


762. 


1009. 






1-3 


0 . 


255. 


511. 


760. 


1007. 


2 


0.061 


4-6, 
NOS PD 2 


3. 


259. 


514. 


762. 


1007. 






1-3 


0 . 


257. 


511. 


760. 


1007. 


3 


0.041 


4-6, 
NOS PD 3 


0 . 


257. 


510. 


759. 


1004. 






1-3 


0 . 


256. 


510. 


760. 


1006. 


4 


0.024 


4-6, 
NOSPD 4 


1 . 


259. 


513. 


762. 


1007. 






1-4 


0 . 


255. 


511. 


760. 


1008. 


5 


0.015 


5-6, 
NOSPD 5 


-2. 


255. 


510. 


759. 


1004. 
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PERFORMANCE EVALUATION OF TCP RADIAL TIJPPINE 
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DUCFD TC STANDAP’ AIP IN ACCORDANCE WITH NASA METHOD 
TAL INLET PRESS. = 1*. 7 PSJA, TOT At. INLET TF M PERATHRI 
518.7 DEG. R., GAW4 = 1.4, C° = 0. ?4 PTU/L BM-DEG. F . 
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APPENDIX C 



PROGRAM NOS PD 



Program NOSPD calculates the referred moments and referred flow 
rates for the stopped-rotor tests and is in essence a much simplified 
version of program RADIAL. The computational techniques used are the 
same in both programs and will not be repeated. Program NOSPD is listed 
in Table Cl, and the input and output date are given in Table C2. Table 
B3 lists the torque calibration data used for each run. The items of 
input data defined in the preceding section apply to this program, 
except as noted below. 

Card No . FORTRAN Description 

1 NRUNS Number of runs (clearances) to be processed 



2 



CLNC 



3 



TCD 



4 



NPTS 



5 



V4, V5, 
TCJ, PATM 
PSP, DPVC 
TQ, TRM 
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TABLE Cl 

LISTING OF PROGRAM NOSPD 
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RESULTS OF ZERO-SPEEO TESTS 
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